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I. INTRODUCTION 



The spin s = \ Heisenberg quantum spin chain describes the magnetic properties of 



The spin dynamics 
From a theoretical 



quasi-one-dimensional antiferromagnetic compounds like KCuF 3 
is experimentally investigated using inelastic neutron scattering 
standpoint, the quantity of interest is the dynamic structure function (DSF) S of two local 
spin operators. This is because the magnetic scattering cross section per magnetic site is 
directly proportional to it lal. 

The Heisenberg model has been studied quite intensively |4|, |5(, and because of the pres- 
ence of the Uq(sl 2 ) symmetry f|, a number of exact results are available. Static properties 
need only the Yang-Baxter relation j^j], whereas dynamic correlation functions require the 
additional notion of vertex operators and exploit bosonization methods l7f . One thus obtains 

n 

compact expressions for form factors [6J. 

Regarding the dynamic structure function, the focus has so far been on S2, the two-spinon 
contribution to the total S. First there has been the Anderson (semi-classical) spin-wave 
theory j^J, an approach based on an expansion in powers of 1/s and hence, exact only in 
the classical limit s = 00. It can describe with some satisfaction compounds with higher 
spins |9i|, but fails in the quantum limit s — |. For this latter system, the Miiller ansatz 
las been proposed, built mainly from finite-chain calculations and symmetry considerations 
3. It is an approximate express™ for S 2 that accounts satisfactory tot tnany aspects of 
the phenomenology |2j. More recently, an exact expression for S2 has been obtained, making 
extensive use of the U^ski) symmetry an d a comparison with the Miiller ansatz shows 
that it gives a better account of the data |l2j, [l3| . 

Beyond the two-spinon DSF S2 are of course all the other S , 2 P >2 contributions. The first 
one to look at is the four-spinon dynamic structure function S4, and the purpose of the 
oresent article is to describe its behavior. An exact expression for S4 has been derived in 
1^ | . see also We intend to use that expression to describe the behavior of S4 as a 

function of the neutron energy transfer uj and neutron momentum transfer k. 

We must mention that a preliminary investigation into the behavior of S4 has already 



been initiated in 



16| . There we have used quadratures to perform the integrations involved 



in the expression of S4. But because of slow convergence of the algorithms we wrote, we 
could describe the behavior of 5*4 only as a function of k, and only for relatively small values 



2 



of oj. In the present work, the integrations are performed using Monte Carlo methods, and 
this makes it possible not only to study 5*4 as a function of k for a wider range of values 
of u>, but obtain a description of S4 as a function of u for a wide range of values of k as 
well. It is however important to note that, for the same values of ui, the behavior of S4 as 
a function of k we obtain here by Monte Carlo methods is similar and consistent with the 
„ M we obtained in H^g^at^. 

Also, we systematically carry a comparison of each result we obtain regarding S4 to a 
corresponding one regarding S^. The reason is that there is good familiarity in the literature 
with the two-spinon DSF and so, such a comparison allows faster acquaintance with the four- 
spinon contribution. To have the discussion as clear as possible, we scale both S4 and S 2 to 
one. All results concerning 5*2 are already known only those concerning S4 are new. 

This article is organized as follows. After these introductory remarks, we describe in the 
next section the Heisenberg model and give the definition of the dynamic structure function. 
We write its decomposition in 2p-spinon contributions and give the expressions of S2 and S4. 
In section three, we first determine the four-spinon continuum, the region in the (k, u) plane 
outside which S4 is identically zero. We will see that it extends beyond the spin-wave des 
Cloizeaux and Pearson boundaries. To the best of our knowledge, this result is a first direct 
and explicit exact theoretical confirmation that the total 5* for the infinite chain tails outside 
the spin-wave continuum. Next in this section is a description of the behavior of S4 as a 
function of u for fixed values of k followed by a comparison with the corresponding behavior 
of Si- It is seen that there is consistency with the four-spinon continuum and that the overall 
shape (not the details) of S4 is similar to the one of S2, each in its own respective continuum. 
Last is carried a description of the behavior of S4 as a function of k for fixed values of to 
and a comparison with the corresponding one of S2. Here too similarity between the two 
overall shapes is found as well as consistency with the four-spinon continuum. Section four 
includes concluding remarks and indicates few directions in which this work can be carried 
forward. 

n 

This paper is a continuation of the work [14]. We use the same notation except for a 
slight modification in f)2. 18j) where here we introduce the function h instead of the function 
/ with the relation / = exp (—h). 
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II. FOUR-SPINON DYNAMIC STRUCTURE FUNCTION 



The antiferromagnetic spin-| XXX Heisenberg chain is defined as the isotropic limit of 
the XX Z anisotropic Heisenberg Hamiltonian: 

1 oo 

H = "2 E A« + i) • (2-1) 

n=— oo 

A = (q + g _1 )/2 is the anisotropy parameter and the isotropic antiferromagnetic limit is 
obtained as A — » — l - , or equivalently g — > — 1 _ . Here cr^** are the usual Pauli matrices 
acting at the site n of the chain. The exact diagonalization of this Hamiltonian is performed 
directly in the thermodynamic limit. This is necessary if we want to exploit the U^s^) 
quantum group symmetry present in the model One consequence is the appearance of 
two vacuum states |0),, i — 0,1 due to two different boundary conditions on the infinite 
chain. The Hilbert space T consists of n-spinon energy eigenstates £ n )ei,...,e»;i such 

that: 

n 

mi,---,Qe 1 ,...,e n ;i = Y,<ZM>--->Oeu..,en-,i, (2-2) 

3=1 

where e(^) is the energy of spinon j and £ ■ is a spectral parameter living on the unit circle. 
In the above relation, ej = ±1. The translation operator T which shifts the spin chain by 
one site acts on the energy eigenstates in the following manner: 

n 

r|^,...,Oe 1 ,.. ) e n ;i = n T ^)l^'---'Oe ll ... ) e n ;l- i , (2-3) 

i=l 

where t(£ •) = e~ ip ^^ and is the lattice momentum of spinon j. The exact expressions 
of the spinon energ y a nd lattice momentum in terms of the spectral parameter are known 
in the literature lall^. \v\ . We are interested in their XXX limit and it is given below in 
eq (j2.12j) . The completeness relation in T reads: 

I 

i=0,l n>0 {e j =±l} j=1 , n "' " 3=1 *~ 

The two-point dynamic structure function is defined as the Fourier transform of the zero- 
temperature vacuum-to-vacuum two-point function. The transverse DSF is therefore defined 
by: 

/oo 
dt ^ t+km) ,<0|a+(tK-(0)|0>, , (2.5) 
-oo 



1 r n rl£ 

EE E ^ tU^\^"^^^^^---u- (2-4) 



where u and k are the neutron energy and momentum transfer respectively, and denotes 
[o x ± ia y )/2. The DSF satisfies the following relations: 

S^iu, k) = S^'ioo, -k) = S^ + -{u, k + 2vr) , (2.6) 

expressing reflection symmetry and periodicity. Inserting the completeness relation (|2.4|) 
and using the Heisenberg relation: 

a x m v ' z (t) = exp (iHt) T' m a X Q y ' z (0) T m exp {-iHt) , (2.7) 

we can write the transverse DSF as the sum of n-spinon contributions: 

S*' + -(uj 1 k)= Stl + -(uJ,k), (2.8) 

n even 

where the n-spinon DSF S n is given by: 

' m& ei,...,e„ J j=l S J \ j=l , 

x Xj£U(£„, -lU , (2-9) 

a relation in which X 1 denotes the form factor: 

In relation (|2.9|) . i + m is to be read modulo 2. Note that each S n satisfies the symmetry 
relations (|2.6jl . 

An exact expression for the form factor X* is known jfJQl- To arrive at the result, one 
has to exploit extensively the infinite-dimensional representation of U g (sl2) and bosonize the 
relevant vertex operators in order to be able to manipulate in a systematic way traces of 
these operators which ultimately yield the correlation functions. Using this form factor, it 



is possible to give an exact expression for the n-spinon DSF in t 



determine its isotropic limit 



obtained via the replacement 



ie anisotropic case [l^ and 
U 



t = ie~ 2%ep ; q = -e~ £ , e - 0+ , (2.11) 

where p becomes the spectral parameter suited for this limit. The expressions of the energy 
e and momentum p in terms of p then read: 

7T 

e (p) = n r = — 7T sin p ; cot p = sinh(27rp) ; — n < p < . (2.12) 

cosh(27rp) 



It turns out that the transverse two-spinon DSF S? does not involve a contour integration, 

3een derive 

1 e~ I{ - p) 

S% (u,k - 7r) = - = 8 (fa; - u 2 i) Q{u>2u - oj) , (2.13) 



see f!2.9|) . Its exact expression has been derived in |ll|. It reads: 



Sod 



where 6 is the Heaviside step function and the function is given by: 

dt cosh(2t) cos(4pt) - 1 t 



Hp)= / \ . , ; . — e\ (2.14) 

; ' t sinh(2t) cosh(t) v ; 



o 



u>2 U (i) is the upper (lower) bound of the two-spinon excitation energies called the des 
Cloizeaux and Pearson (dCP) 10, 11] upper (lower) bound or limit. They read: 

k 

oJiu = 2tt sin — ] uj 2 i = it | sin k\ . (2.15) 

The quantity p is related to u and k by the relation: 











u 2 - 





cosh7rp=W^ f, (2.16) 



which is obtained using eq (12.121) and the energy-momentum conservation laws. The prop- 

nn n 

erties of 5*2 have been discussed in |12L [13J where a comparison with the Miiller ansatz [10J 
is carried. 

The four-spinon DSF S4 involves only one contour integration and its expression is given 
in For < k < n it reads: 

rfp 3 / d Pi F( Pl ,...,p 4 ). (2.17) 

For other values of k, it extends by symmetry using ()2.fij) . C4 is a numerical constant 
irrelevant for the present work since we will scale 5*4 to unity, and the integrand F is given 
by: 

F( Pu ...,p A )- Y -,P4)] EU\9e( Pl ,..,pd\\ (21g) 

The different quantities involved in this expression are as follows. 

W = uj + 7r (sin p 3 + sin p 4 ) ; 
W u = 2vr|sin (AT/2) | , K = k + p 3 + p 4 ; 
cot Pj = sinh(27rp J ) , — 7r < < . (2-19) 
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The function h is given by: 



% 1; ..., P4 )= Yl J (%')' ( 2 - 2 °) 

l<7</<4 



where p„/ = p ■ — p •/ and the function gi reads: 

4 

» = (-l) m (2 7 r) 4 ^cosh(2vrp J ) 

3=1 

X m= ty_ £) IL^ ^ sinh(vrp J ,) 11 F (m + 1 + ip^) ' 

where is the Heaviside step function. In (|2.18|) . the sum ^(p lP2 ) * s over the ^ wo P a i rs 
(pi,P2) and (p2,Pi) solutions of the energy-momentum conservation laws: 

= — 7r(sin pi + sin p 2 ) ; K = — pi — p 2 . (2.22) 

They read: 

( Pl ,p 2 ) = (-iT/2 + arccos (W/ [2tt sin (i^/2)]) , -AT/2 - arccos (W/ [2tt sin (AT/2)])) . 

(2.23) 

Note that the solution in (|2.23|) is allowed as long as W\ < W < W u where W u is given in 
flTiTH) and: 

Wi = ?r| sin K\ . (2.24) 

The (analytic) behavior of the function F in (|2.18|) is discussed in Q|. It is shown that 
the series g% is convergent. It is also shown that gi stays finite when two p^s or more get 
equal. Since the function exp (—h) goes to zero in these regions the integrand F of S4 
is regular there. Furthermore, it is shown that F is exponentially convergent when one of 
the p/s goes to infinity, which means the two integrals over p 3 and p 4 in (j2.17|) do not yield 
infinities. All these analytic results help secure safe numerical manipulations. 



III. BEHAVIOR OF EXACT FOUR-SPINON DSF 

From now on, we restrict ourselves to the interval < k < it. All forthcoming results can 
be carried to the other intervals of k using the symmetry relations (|2.fij) . Also, we scale both 
S4 and S2 to appropriate units in order to display conveniently their respective behaviors. 



7 



A. Four-spinon continuum 



The first feature we discuss is the 'four-spinon continuum', by analogy with the two- 
spinon (or the spin-wave) continuum. It is the extent of the region in the (k, a»)-plane 
outside which S4 is identically zero. Remember that from (|2.13|) . S2 is confined to the 
region uj2i(k) < uj < U2 U {k), where 0J2i, u {k) are the dCP boundaries given in (j2.15|) . From 
the condition W\ < W < W u mentioned after ()2.23|) . we deduce that in order for S4 to be 
nonzero identically, we must have 0Ju(k) < uo < U4 U (k), where: 

^4i{k) = 37rsin(fc/3) for 0<k<n/2] 

&4i(k) = 37r sin(fc/3 + 27r/3) for 7r/2<k<7i; 

uj Au (k) = 4vrcos(A;/4) for < k < n . (3.1) 

We see that u>u(k) and uj4 U (k) are a sort of four-spinon dCP boundaries for S4. The two 
and four-spinon continua are drawn in figure 1 below. We immediately notice that the 
four-spinon continuum is not restricted to the region between the two-spinon dCP branches, 
which means, a fortiori, that the full 5* is also not confined to the spin-wave continuum. 
This is a direct and explicit theoretical confirmation of the 'tail' of the dynamic structure 
function observed outside the spin-wave continuum in finite-chain numerical calculations 
and the phenomenology Q]. Due to the six and higher spinon contributions, it is actually 
legitimate to expect the full DSF to tail even further outside the four-spinon continuum, 
with arguably much smaller values. For example, in the interval < k/ir < 0.51741, there 
is a narrow region between U2 U and inside which both S2 and S4 are identically zero 
whereas the total S may have (very small) nonzero values, something that could eventually 
be checked in finite-chain calculations. However, it is not possible to estimate exactly the 
continua corresponding to the S n> 4 without manageable explicit formulae for these. But 
what is already clear from our study is the fact that indeed, the spin-wave continuum is not 
restrictive to the total dynamic structure function. 

There are further general features we can also read from figure 1 without having recourse 
to detailed calculations. For example, we see that inside the interval < k/71 < 0.51741, 
the four-spinon continuum lies entirely above the two-spinon continuum. This means that 
for this interval, S2 may be accepted as a good approximation for the total S between the 
spin-wave boundaries and 64 between uj^ and u>^ u . For 0.51741 < k/n < 1 however and as k 



S 



FIG. 1: Four-spinon (full) and two-spinon (dashed) continua. 



increases, there is increasing overlap between the two continua so that inside the spin-wave 
continuum, we may expect S to divert in perhaps a non- negligible way from 5*2. 

B. Behavior as a function of the energy transfer 

Next we describe the behavior of S4 as a function of uo for fixed values of k. Figure 2 
displays the shapes of S4 for respectively k/n = 1/4, 1/2, 3/4 and 1. As an illustration, take 
for example the case k/n = 1/4. We see on figure 2 that S4 is zero until uj/tt ~ 0.68, which 
corresponds to the beginning of the four-spinon continuum u>4i(ir/ '4) /n = 3sin(7r/12) = 
0.77646. We obtain a slightly smaller value when reading directly from figure 2 because 
of fitting. Note that it is important to check each time consistency with the four-spinon 
continuum. This is because the lower and upper branches u>4i ;U (k) are not imposed explicitly 
in the integration algorithm of S4 unlike the case of £2 where the corresponding expression 
(|2.13j) incorporates explicitly the two-spinon dCP boundaries 0J2i, u (k). For S4, only the 
conditions W\ < W < W u noted just before (|2.24j) are imposed. 

Starting from uo/Att ~ 0.17, S4 jumps sharply from zero to a maximum at uj/Att ~ 0.21 
(read from figure 2). Then it decreases with two apparent local minima until it becomes 
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FIG. 2: (Scaled) S4 as a function of to for fixed k. 



negligible at roughly u/An ~ 0.6. The upper branch of the four-spinon continuum at k = ir/4 
is u; 4m (7r/4) /47r = cos(7r/16) = 0.9808. Hence, S4 is practically negligible (but not identically 
zero) for values of uj/Att between about 0.6 and 0.9808. The description of S4 as a function 
of uj for the other values of k can be carried along the same lines and each time, consistency 
with the four-spinon continuum is checked. The shapes are all similar to one another: steep 
increase from zero to a maximum value, followed by a 'wiggled' slower decrease to zero. 

Furthermore, it is interesting to notice that for all values of k, the overall shape of 
S4, not the detail, is roughly similar to the one of S 2 , represented as a function of uj for 
the same values of k in figure 3. For instance, for /c/7r = 1/4, S2 starts very sharply 
at about u/2n ~ 0.35, which corresponds to the start of the two-spinon continuum at 
^2«( 7r /4)/27r = 0.35355. It reaches a maximum before decreasing more slowly towards zero. 



10 



FIG. 3: (Scaled) S2 as a function of to for fixed k. 



The only difference worth mentioning is the decrease of S4 after its absolute maximum which 
presents local minima and maxima whereas the decrease of S2 is always 'smooth'. This may 
originate from the 'richer' structure of the expression of S4 with respect to that of S2. 

C. Behavior as a function of the momentum transfer 

Last we describe the behavior of S4 as a function of k for fixed values of uo. In figure 
4 are plotted the graphs of £4 in terms of k for uj/tt = 1/4, 1/2, 3/4, 1 respectively. Let 
us describe for example the case uj/tt = 1/4. We see on figure 4 that £4 is zero until we 
reach the value k/n ~ 0.88. On the other hand, the four-spinon continuum lies outside the 
interval 0.07967 < k/n < 0.92033. In the region < k/n < 0.0796 7, figure 4 shows no 
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discernible finite values for S4, only a very thin 'trace' that would be more visible with a 
better resolution. This means that S4 is negligible for those small values of k. For larger 
values of cu though, S4 picks up clear finite values in the interval < k < 3 arcsin(o;/37r); 
see the other graphs on figure 4. Those values get larger as uj increases. 



FIG. 4: (Scaled) £4 as a function of k for fixed uj. 

Returning to the case uj/ti = 1/4, we see that £4 starts from zero at k/n ~ 0.88 (slightly 
smaller than the exact value 0.92033 because of fitting), rises sharply to a maximum and 
then decreases. The behavior of S4 for the other values of uj is also consistent with the 
four-spinon continuum. Take for example the case uj/tt = 1/2. The four-spinon continuum 
indicates that is identically zero for 0.1599 < k/n < 0.84008. We see indeed small values 
for ^4 from k/n = up to a little before 0.2, and rising again from zero a little after 
kjix = 0.8 to a local maximum, then to an absolute maximum before decreasing. 
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FIG. 5: (Scaled) 5*2 as a function of k for fixed u>. 



Here too S2 has a similar overall behavior within its own (two-spinon) continuum. Figure 
5 displays the behavior of S2 with respect to k for the same fixed values of to. But as u 
increases, we notice the richer structure of 5*4 with respect to the corresponding one for 
S 2 - This is due to the more involved expression of S4. In fact, for larger values of u>, the 
integrand F in (|2.17|) is nonvanishing in larger and larger areas in the (p3,Pa) plane. For 
illustration, compare the behavior of F shown in figure 6 for (k, u) = (0.57T, 2tc) with the one 
shown in figure 7 for (Lw) = (0.57T, 3n). In this regard, we recall that the quadrature-based 
algorithms written in [l6| were based on this observation and took advantage of the fact 
that for fairly small values of u, the integrand F was negligible in large areas of the (p3,Pa) 
plane. This is the reason why those algorithms could not be carried to larger values of ui, 
or be able to describe efficiently a behavior as a function of the neutron energy u itself. As 
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FIG. 6: The integrand F in ( I2.17|) as a function of (p3,P4,) for (k,u) = (0.52-7T, 2ir). 



FIG. 7: The integrand F in ( I2.17|) as a function of (p3,p^) for (k,u>) = (0.57r,37r). 

already mentioned, we have used here Monte Carlo techniques and we can do better, but 
consistency with |l6( is realized, i.e., we have the same behavior of S4 as a function of k for 
the same small values of uo. 

IV. CONCLUSION 

In this work, we have described the behavior of the exact four-spinon dynamic structure 
function S4 in the antiferromagnetic isotropic Heisenberg quantum spin chain at zero tem- 
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perature as a function of the neutron energy u and momentum transfer k. We have also 
determined the four spinon continuum, the region outside which S4 is identically zero. The 
discussion was carried in the form of a comparison with the corresponding behavior of the 
exact two-spinon dynamic structure function S2, already known in the literature. Figures 8 
and 9 summarize these two behaviors where S4 (8) and S 2 (9) are drawn in the (k, ui) plane. 
Recall that we have scaled them down to one. 



FIG. 8: (Scaled) S4 as a function of k and u. 



FIG. 9: (Scaled) S2 as a function of k and u. 
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There are four directions in which one may wish to carry forward with this work. The 
first is the anisotropic case. The model is exactly soluble and we do have generic expressions 
for S n in the form of contour integrals in the spectral parameters complex planes H- The 
difficulty here is that the integrands involve much more complicated functions which are 
already present in S2, and one should expect intricate complexities in this more general 
case. 

The second direction in which one may want to push forward is the situation where there 
is an external magnetic field. There are finite-chain calculations in this regard, and more 
recently the works [19]. But one has to remember that the model in not exactly solvable in 
this case. So one may want to try small perturbations around the zero-field limit solution. 
The third direction is the finite-temperature case. Here too there are finite-chain results and 
it is interesting to see those effects on S2 and S4. The fourth direction is to look into the 
situation of a spin-one chain. The model is still exactly solvable and exploiting the quantum 
group symmetry, compact expressions for the form factors are available [2^ , 3, . 



[1] K. Hirakawa and Y. Kurogi, Prog. Theor. Phys. S46 (1970) 147. 

[2] D.A. Tennant, R.A. Cowley, S.E. Nagler and A.M. Tsvelik, Phys. Rev. B52 (1995) 13368; 
D.A. Tennant, S.E. Nagler, D. Weltz, G. Shirane and K. Yamada, Phys. Rev. B52 (1995) 
13381; D.A. Tennant, T.G. Perring, R.A. Cowley and S.E. Nagler, Phys. Rev. Lett. 70 (1993) 
4003; S.E. Nagler, D.A. Tennant, R.A. Cowley, T.G. Perring and S.K. Satija, Phys. Rev. B44 
(1991) 12361. 

[3] G.L. Squires, 'Introduction to the theory of thermal neutron scattering', Cambridge University 
Press, 1996; S.W. Lovesey, 'Theory of neutron scattering from condensed matter', Clarendon, 
Oxford, 1987. 

[4] Some (not all) of the literature related to the Heisenberg model comprises W. Heisenberg, Z. 
Phys. 49 (1928) 619; H. Bethe, Z. Phys. 71 (1931) 205; L. Hulthen, Arkiv. Mat. Astron. Fysik 
A 1126 (1938) 1; E.H. Lieb and D.C. Mattis, J. Math. Phys. 3 (1962) 749; J. des Cloizeaux 
and J.J. Pearson, Phys. Rev. 128 (1962) 2131; R.B. Griffiths, Phys. Rev. 133 (1964) A 768; 
C.N. Yang and CP Yang, Phys. Rev. 150 (1966) 321; 150 (1966) 327; 151 (1966) 258; Th. 
Niemeijer, Physica 36 (1967) 377; E. Barouch, B.M. McCoy and D.B. Abraham, Phys. Rev. 



16 



A4 (1971) 2331; M. Gaudin, Phys. Rev. Lett. 26 (1971) 1301; M. Takahashi, Prog. Theor. 
Phys. 46 (1971) 401; L.A. Takhtajan and L.D. Faddeev, Russ. Math. Surveys 34 (1979) 11; 
B.M. McCoy, J.H.H. Perk and R.E. Shrock, Nucl. Phys. B220 (1983) 35; O. Babelon, H.J. de 
Vega and CM. Viallet, Nucl. Phys. B220 (1983) 13; G. Miiller and R.E. Shrock, Phys. Rev. 
B29 (1984) 288; J.M.R. Roldan, B.M. McCoy and J.H.H. Perk, Physica 136A (1986) 255; 
V.E. Korepin, A.G. Izergin and N.M. Bogoliubov, 'The Quantum Inverse Scattering Method 
and Correlation Functions', Cambridge University Press, 1993; F.H.L. Essler, H. Frahm, A.G. 
Izergin and V.E. Korepin, Comm. Math. Phys. 174 (1994) 191; V.E. Korepin, A.G. Izergin, 
F.H.L. Essler and D. Uglov, Phys. Lett. A190 (1994) 182. 

[5] R.J. Baxter, 'Exactly Solved Models in Statistical Mechanics', Academic Press, 1982. 

[6] M. Jimbo and T. Miwa, 1 Algebraic Analysis of Solvable Lattice Models', American Mathemat- 
ical Society, 1994. 

[7] O. Davies, O. Foda, M. Jimbo, T. Miwa and A. Nakayashiki, Comm. Math. Phys. 151 (1993) 
89; LB. Frenkel and N.H. Jing, Proc. Natl. Acad. Sci. 85 (1988) 9373; A. Abada, A.H. 
Bougourzi and MA. El Gradechi, Mod. Phys. Lett. A8 (1993) 715; A.H. Bougourzi, Nucl. 
Phys. B404 (1993) 457; A.H. Bougourzi, ' Bosonization of quantum affine groups and its ap- 
plication to the higher spin Heisenberg model', q-alg/9706015. 

[8] PW. Anderson, Phys. Rev. 86 (1952) 694. 

[9] M.T. Hutchings, G. Shirane, R.J. Birgeneau and S.L. Holt, Phys. Rev. B5, (1972) 1999. 
[10] G. Miiller, H. Thomas, H. Beck and J.C. Bonner, Phys. Rev. B24 (1981) 1429. 
[11] A.H. Bougourzi, M. Couture and M. Kacir, Phys. Rev. B54 (1996) 12669. 
[12] M. Karbach, G. Miiller and A.H. Bougourzi, ' Two-spinon dynamic structure factor of the 

one- dimensional S = 1/2 Heisenberg antiferromagnet' , cond-mat/9606068. 
[13] A.H. Bougourzi, M. Karbach and G. Miiller, 'Exact two-spinon dynamic structure factor of 

the one- dimensional s = 1/2 Heisenberg- Ising antiferromagnef, cond-mat/9712101. 
[14] A. Abada, A.H. Bougourzi and B. Si-lakhal, Nucl. Phys. B497 [FS] (1997) 733. 
[15] A.H. Bougourzi, Mod. Phys. Lett BIO (1996) 1237. 

[16] A. Abada, A.H. Bougourzi, B. Si-Lakhal and S. Seba, ' Four-Spinon Dynamical Correlation 

Function in Isotropic Heisenberg Model', cond-mat/9802271, unpublished. 
[17] L.D. Faddeev and L.A. Takhtajan, J. Soviet. Math. 24 (1984) 241. 

[18] FA. Smirnov, 'Form Factors in Completely Integrable Models of Quantum Field Theory' 

17 



World Scientific, Singapore, 1992. 
[19] M. Karbach, D. Biegel and G. Miiller, ' Quasiparticles governing the zero-temperature dynamics 

of the ID spin- 1/2 Heisenberg antiferromagnet in a magnetic field', cond-mat/0205142; M. 

Karbach and G. Miiller, l Line shape predictions via Bethe ansatz for the one dimensional 

spin- 1/2 Heisenberg antiferromagnet in a magnetic field' , cond-mat/0005174. 
[20] M. Idzumi, Int. J. Mod. Phys. A9 (1994) 4449; 1 Correlation functions of the spin 1 analog of 

the XXZ model', hep-th/9307129. 
[21] A.H. Bougourzi and R.A. Weston, Nucl. Phys. B417 (1994) 439. 

[22] A.H. Bougourzi, l Bosonization of quantum affine groups and its application to higher spin 
Heisenberg model' , q-alg/9706015. 



18 



00=271 



00=3% 




